Laser cooling of a nanomechanical resonator mode to its quantum ground state 
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We show that it is possible to cool a nanomechanical resonator mode to its ground state. The 
proposed technique is based on resonant laser excitation of a phonon sideband of an embedded quan- 
tum dot. The strength of the sideband coupling is determined directly by the difference between the 
electron-phonon couplings of the initial and final states of the quantum dot optical transition. Possi- 
ble applications of the technique we describe include generation of non-classical states of mechanical 
motion. 



In contrast to bulk semiconductors where lattice vi- 
brations form a reservoir with a continuous spectrum, 
nano-scale semiconductor structures support sharp me- 
chanical resonances with very high quality factors and 
frequencies approaching 1 gigahertz [jj. Unlike optical 
resonators, the thermal occupancy of these mechanical 
resonators is well above unity, even when they are cooled 
to sub-Kelvin temperatures. Presence of thermal noise 
hinders the study of coherent quantum dynamics in these 
mechanical systems and limits their applications to 
precision measurements Here, we propose a technique 
that would realize laser cooling of a nano-mechanical res- 
onator mode to its motional ground-state where quantum 
effects will be manifest. The basic idea is to use light scat- 
tering and electron-phonon interactions 4| , such as defor- 
mation potential coupling 5], in an embedded quantum 
dot to manipulate a discrete mode of lattice vibrations. 

The semiconductor beam structure with an embedded 
quantum dot (QD) that we analyze is shown in Fig. 1. An 
infinite length beam supports four phonon branches with- 
out an infrared cutoff: two bending branches (in-plane 
bending and flexural) with quadratic dispersion relations, 
and a torsional and a compression branch with linear dis- 
persion relations @- For a finite beam of length L 
attached to supports via abrupt junctions, these branches 
at long wavelengths will develop a series of sharp reso- 
nances corresponding to the different harmonics ja| . In a 
structure where the thickness of the beam (d) is smaller 
than its width (b) (Fig. 1), the lowest-energy resonance 
corresponds to the fundamental flexural (vertical bend- 
ing) mode (Hujq). This will constitute the resonator mode 
we intend to address and cool to the ground state. For 
flexural modes with uj ~ 10 9 s _1 , quality factors (Q) ex- 
ceeding 20, 000 have been measured [||. Though we focus 
on a bridge geometry (Fig. 1), the analysis of a cantilever 
structure is completely analogous. 

The presence of a zero-dimensional (anharmonic) emit- 
ter that can be resonantly excited by a laser, has a spon- 
taneous emission broadened optical transition, and inter- 
acts with the acoustic modes is necessary for achieving 
ground-state cooling. Our discussion will primarily focus 
on self-assembled InAs QDs .10] which can satisfy these 



requirements and can be modeled as two level systems 
consisting of the empt y Q D state (\g)) and the fundamen- 
tal exciton state (|e))[ll|. We envision an ideal scenario 
where processing leaves only one QD inside the beam at 
a specific location. In assuming a natural line- width lim- 
ited QD we are ignoring the effect of the surface on the 
electronic properties of the dot. This is warranted pro- 
vided the dephasing and/or non-radiative recombination 
rates arising from the proximity to semiconductor/air in- 
terfaces are negligible compared with the spontaneous 
emission rate. 

We observe in Fig. 1 that flexion induces extensions 
and compressions in the structure 0. This longitu- 
dinal strain will modify the energy of the electronic 





QD 




















II 

L 



WL 















III 

L 



FIG. 1: Top: Schematic diagram of a GaAs bridge of length 
L, width b and thickness d (see II below), with an embedded 
lens shaped InAs QD (WL denotes the wetting layer). Bot- 
tom: Vertical cross-sections through the axis of the bridge il- 
lustrating the deformations suffered by a small neighborhood 
of the QD, when the beam is bended slightly in the vertical 
direction (x). I and III correspond respectively to positive 
and negative deflection and II shows the equilibrium configu- 
ration. Each volume clement parallel to the axis undergoes a 
simple extension or compression^. 
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states confined in the QD through deformation poten- 
tial coupling 5]. The resulting interaction between a lo- 
calized carrier and a discrete flexural phonon mode will 
lead to the appearance of sidebands in optical spectra, 
originating from phonon assisted photon absorption or 
emission processes By tuning a laser field into res- 
onance with the first lower energy (red) sideband of a 
QD optical transition, we can ensure that absorption of 
a laser photon will be accompanied by the removal of a 
phonon from the svstem|121 yjfl. 

Before proceeding we note that sideband cooling has 
been used to cool trapped ions to their motional ground 
state 12]. Despite the analogy, the physics of the two 
systems are different. Unlike the case of trapped ions 
[13J and in marked contrast to optomechanical cooling 
schemes 0], the photon recoil plays no role in the laser 
coo ling of a nanomechanical mode using an embedded 
QD[2s|. In the latter system, the (oscillator) strength 
of the sidebands is determined directly by the difference 
between the clcctron-phonon couplings of the initial and 
final states of the optical transition^]. 

The fact that the modes of interest correspond to 
acoustic phonons with typical wavelengths X p ~ L ^> 
b, d (Fig. 1) orders of magnitude larger than the lat- 
tice constant makes thin rod elasticity theory an excel- 
lent approximation to describe the phonons inside the 
beamU0,0]- The dominant electron-phonon interaction 
mechanism in our structure will be deformation potential 
coupling [24| 



H DP = J df 3 [D c p cl (f) - D v p h (f)] V • Q(f) ■ 



(1) 



Here, we have introduced the electron (hole) density 
operator p c \ (p\) and the corresponding deformation 
potential constants D c (D v ). The vector field opera- 
tor u(f) corresponds to the displacement in Lagrangian 
elasticity 0. We neglect the difference between the elas- 
tic properties and deformation potentials of the strained 
InAs QDs and the surrounding GaAs matrix [2^ . 

To proceed we decompose u(f) in terms of a set of 
orthonormal modes. A possible set consists of dis- 
crete modes localized in the beam ({u m (f)}) satisfying 
clamped boundary conditions at the junctions, and a con- 
tinuum of modes localized in the supports ({%(?)}) sat- 
isfying free boundary conditions. In this representation 
the beam modes only interact with support modes [lfij 
and for long wavelengths (X p ) these interactions are weak, 
making the {u m (f)} correspond to the resonances 0. As 
we are just interested in the dynamics of the fundamental 
mode (described by uo(f), 60), we single it out and diag- 
onalize all the terms in the phonon Hamiltonian which 
do not involve bo, by transforming the rest of the modes 
to a new continuum {u q (f)}. Thus we obtain a represen- 
tation consisting of the discrete resonator mode weakly 
coupled to a reservoir of "background" modes {b q }. 



We subsequently restrict JIJ to the two level model for 
the QD and integrate over f to obtain a coupling term 
between the QD and the resonator mode of the form: 
hojoi]\e) (e\(b Q + b ), where (Do is the bare frequency asso- 
ciated to the discrete resonator mode. Within thin rod 
elasticity, V • u(f) for a bending mode is well approxi- 
mated by (2cr - l)x^r, where X(z) (Fig. 1) is the de- 
flection of the center-of-mass of the cross-section of the 
beam, a is the Poisson ratio Q and the coordinate ori- 
gin is at the midpoint. We use this relation, that A p is 
much larger than the typical QD size 11], and that the 
differences between the average positions of the hole, the 
electron and the center of mass of the QD can be ne- 
glected, to obtain: 

2 _ 3V6(l-2a) 2 (D c - D v ) 2 x%L xf (k z D ) 

^ — „. o ,3/2 J4I. r / n2 • \ l > 



^hpcip (1 + er) 



(koL/nY 



To the extent that one neglects the difference between the 
bare frequency ujq and the actual resonance ujq, tj 2 cor- 
responds to the Huang-Rhys parameter of the resonator 
mode. Here p and ct are the density and transverse speed 
of sound for the material of the beam, and and zjj are 
the QD coordinates (Fig. 1). The mode function normal- 
ized to the length and the wave-vector of the resonator 
mode are denoted by Xq (koz) and fco, respectively. The 
first factor in Eq. (2) is the square of a characteristic 
length of the material (2 nm for GaAs). The second 
factor gives the dependence with the dimensions of the 
structure and the placing of the QD in the vertical direc- 
tion. The last factor will depend on whether we have a 
cantilever or a bridge, and for the latter is of order unity 
in a neighborhood of the midpoint. Typical parameters 
(see below) yield 77 ~ 0.06. 

The Hamiltonian describing our system includes terms 
corresponding to the interactions between the QD, the 
resonator mode, the background phonon modes, and the 
electromagnetic field. To study laser cooling we per- 
form a canonical transformation that eliminates the QD- 
resonator mode coupling termy. After additional co- 
ordinate and momentum shifts applied to the phonon 
modes, we obtain the transformed Hamiltonian: 



H = H - hS- 



h.c. 



, (3) 



where H = h[Cj b f b + J2 q u q b\b q + J2 k i^k ~ wl) a\a k \ 
and the QD operators are expressed in terms of Pauli ma- 
trices (i.e. |e)(e| = l/2(l+<7z)). We have also defined the 

polaron operator B = e r, ( b °~ b °\ the detuning of the laser 
from the fundamental exciton line S and its frequency cul 
and Rabi frequency Q. The background phonon modes 
are characterized by annihilation operators b q and their 
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couplings to the QD (resonator mode) X q ((g). The radi- 
ation field is characterized by annihilation operators ak 
and their couplings to the QD g^. 

The role of r\ in our system is equivalent to the 
one played by the Lamb-Dickc (LD) parameter in the 
trapped-ion svstem[Hll7t. Even though the underlying 
physics is different, we will refer to the LD regime as the 
one in which rj 2 s = r) 2 ((blb Q ) + 1) 1. We focus on 
this limit since it is the one of interest for ground state 
cooling and is well satisfied for realistic structures pro- 
vided the ambient temperature is low enough. In this 
regime, both cooling and heating can take place via two 
paths, where the modification of the motional state of 
the phonon mode takes place in the laser absorption and 
spontaneous emission process, respectively (Fig. 2). In 
contrast to the trapped ion system, these two paths are 
indistinguishable for the QD-beam structure and exhibit 
quantum interference. 

In this LD regime it is possible to derive a rate equa- 
tion describing the cooling dynamics of the resonator 
mode. To this effect we first eliminate the radiation field 
and subsequently the "background" phonons to obtain a 
master equation for the QD-resonator system. The first 
step, that involves the standard Born-Markov approx- 
imation, leads to a Liouvillian for the "QD-resonator- 
background" system that - aside from terms involving 
the "background" - closely resembles the analogous re- 
sult for a trapped- ion[l^|. The only difference -aside from 
a trivial interchange of the normal coordinate and mo- 
mentum of the resonator mode- lies in the dissipative 
part: T/2(2a-BpB^a + — pa + er_ — er + (7_ p) , where p is the 
density matrix. The crucial point is that B does not de- 
pend on the electromagnetic modes, in marked contrast 
to the momentum shift operator that plays the analogous 
role for a trapped-ion. This gives rise to the interference 
effects already discussed. 

The background will exhibit sharp resonances corre- 
sponding to the other beam modes (u m (f)\ m ^o). It can 
be shown that provided these resonances are well detuned 
from ujq and 2ojq (see H3)) the effect of the couplings X q 
in Hamiltonian J3J on the cooling dynamics of the res- 
onator mode is higher order in the small parameters rj 2 
and d/L (Ref. ^H). Hence we neglect these couplings in 
our derivation of a master equation for the QD-resonator 
system. The weakness of the interaction between the 
resonator and the supports allows us to apply the 
Born-Markov approximation to the remaining couplings 
{Cq}- This yields three distinct contributions: a Hamilto- 
nian part that shifts Qq to u>q, a "counter-rotating" part 
that to lowest order in 1/Q will not contribute to the 
rate equation, and a "rotating wave" part describing the 
damping of the resonator [fjj. 

We proceed to apply the LD approximation following 
Ref. |17| . This involves expanding all the polaron ex- 
ponentials B up to second order in rj and adiabatically 
eliminating the QD. The damping of the resonator can 
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FIG. 2: Energy level diagram of the QD coupled to the res- 
onator mode for perturbative Rabi frequency. \g,n) (|e, n)) 
denote the state where the QD is in the ground (exciton) 
state and the resonator mode has n phonons. The cooling 
(heating) cycle consists of the excitation of the QD and sub- 
sequent spontaneous photon emission, and is denoted by solid 
(dashed) lines. The blue (red) path corresponds to the case 
where phonon number is changed during the laser absorption 
(spontaneous emission) . If the laser is tuned between the red 
sideband and the main line these two paths interfere con- 
structively for cooling and destructively for heating (note the 
relative phase of n between vs and «l). 



be included using an "LD-large Q" approach in which 
1/Q is treated on the same footing as rj 2 . The resulting 
rate equation describing the populations P n of the energy 
levels of the resonator mode is 



Pn = 



CJQ- 



Q 



rfA- 



UJQ 



»(^o) + l 

Q 



[nP„_i - (n + l)P„] 

[{n + 1) P„+i - nP n ] , (4) 



with n(uJo) = l/[e Wo / fe£iT — 1]. The first term corresponds 
to heating and the second one to cooling. There are 
contributions proportional to 1/Q describing the damp- 
ing arising from coupling to the supports and contribu- 
tions which involve the scattering of laser light given by 
ri 2 A± = 2r, 2 Re{J2„A ± 0lJ (a„) ss }, with 



A ± = ^C+(f W -TC)iTiu I-~Ly 1 -(% R»-rc). 

(5) 

We have used for the state space of the QD a basis 
consisting of the Pauli matrices (a v with v > 0) and 
the identity (er ), and defined: R y v = |Tr QD (a^ayCT^), 
C^u = |Tr QD (a f _ l a-a u a+) and L M „ = |Tr QD (cr M £ QD ov). 

Here £ QD and (<t i/ )ss ar e the Liouvillian and steady state 
expectation values corresponding to the optical Bloch 
equations for the QD0|. We focus on S < for which 
there is a net cooling rate r] 2 W = r} 2 (A^ — A + ) > 0. 

As we start from thermal equilibrium, initially the 
number of phonons is given by rn — n(tuo). When steady 
state is reached, the final phonon number n/ will be given 
by nf = (m + rfQA + /uj )/(l + i] 2 QW/uj )- Therefore 
there is a threshold value of = min{A + /H^}|n il 5 below 
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which laser cooling does not work (Fig. 3), and rj 2 Q >> 1 
is required to obtain appreciable cooling. 

There are two contributions ton/: the first one is pro- 
portional to m and the second one proportional to A+. 
This leads to two regimes depending on which contri- 
bution dominates the behavior of nf (Fig. 3). In the 
first "large n^-small rj 2 Q" regime the dominant heating 
mechanism is the coupling to the supports. In the second 
"small n^-large i] 2 Q" regime heating is dominated by the 
scattering of laser light and the optimal value of rif (hf) 
is given by hf ~ (r/4w ) 2 (Fig. 3). We find that in both 
regimes T /ujo < 1 is a necessary condition for ground 
state cooling and hf < Ui/rfQ + (r/4wo) 2 - Typical pa- 
rameters we envisage are: L = 950nm, d = 30nm (which 
correspond to uiq w 1.2 x 10 9 s _1 ), b = 85nm, xd = d/4, 
z D = 0, Q = 3 x 10 4 and V = 3 x lO^"^- These lead 
to hf = 0.1 for a sample held at T = 0.1K (Fig. 3). 

We note that it is possible to realize an all optical 
measurement of the final temperature. The basic idea 
is to take two lasers, a first one for cooling with the de- 
sired parameters, and a second one for measuring with 
zero detuning (frequency u>m) and perturbative Rabi fre- 
quency (e.g. Q,m = r/10). While this second "probe" 
laser is always on, the first "cooling" laser interacts with 
the QD only during time windows of duration t c (cool- 
ing) that alternate with measurement windows of dura- 
tion t m -C t c . As for realistic parameters optimal cool- 
ing occurs at Q ~ luq, when analyzing the cooling win- 
dows the "probe" laser can be neglected. The experi- 
mentally measured quantities are the photo-counts iV_ 
(blue photons) and N+ (red photons) around frequen- 
cies u>m + and ujm — respectively, during the time- 
intervals when the "cooling" laser is off. The cooling and 
heating rates satisfy A T (S) = A±(— S). Hence if we as- 
sume that there are no background photons, it follows 
that in the above scenario the occupancy will be given 
by(btb )=N_/(N + -N_). 

We have focused on laser cooling the fundamental res- 
onator mode to its ground state. However we want to 
stress that one of the main contributions of this letter 
is to establish a close analogy between this mesoscopic 
condensed-matter system and a single trapped ion in the 
LD regime 0>0|. In fact, the analogy extends to cavity- 
QED; i.e. a single atom in a high-Q optical or microwave 
cavity IfSI. This opens up possibilities that go far beyond 
cooling |20j and may not be within reach of other tech- 
niques for manipulating a nano-resonator. In particular 
the technique presented here would enable quantum state 
engineering of non-classical states of motion 0, Il3l along 
the lines proposed for the trapped-ion system [21j. Ex- 
amples are the generation and detection of a Fock state 
and squeezed states of motion. 

One could also envision extending this method to si- 
multaneous cooling of several modes. If successful, this 
could be used to generate entangled states of mechan- 
ical motion. Our results could apply to other systems 
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FIG. 3: Optimal final number of phonons hf v.s. initial num- 
ber m (solid line) for F/ojo = 1/4, r) 2 Q = 100. The optimum 
results from searching for the laser parameters (5 op t, f2 op t) 
that minimize nf. There are three regimes (dotted lines): I) 
for m < (g^j) 2 laser addressing of the structure is not useful 
(i.e. rif > rii), II) for small m (rii <§C ifQij^) 2 ) the op- 
timum is constant and <5 op t ~ Wo, ^o P t "C Wq; HI) f° r large 
rii (rii 2> V 2 Q(^^) 2 ) the optimum is proportional to m with 
the ratio hf/rii bounded by l/r] 2 Q and 5 op t ~ <^o, ^o P t ~ Wo. 
The rough approximation ni/r] 2 Q + (r/4ojo) 2 is plotted for 
comparison (dashed line). 

with sharp vibrational resonances and to different t ype s 
of zero-dimensional emitters such as localized defects |22| . 
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